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^ ■ Large-scale magnetic fields can affect scalar cosmological perturbations whose evolution is 

described in the conformally Newtonian gauge and within the tight coupling approximation. 
>. ■ The magnetized curvature perturbations present after matter radiation equality (and prior 

j_j ■ to decoupling) are computed in terms of an appropriate transfer matrix allowing a general 
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estimate of the Sachs- Wolfe plateau. From the observation that CMB initial conditions 
should be (predominantly) adiabatic, the contribution of the magnetic field intensity can be 
constrained. 
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Large-scale magnetic fields are observed at a /i G level in galaxies, clusters and in some 
superclusters [1]. Compressional amplification (taking place during the gravitational collapse 
of the protogalaxy) allows to connect the observed magnetic field to a protogalactic field, 
present prior to gravitational collapse, of typical strength of 0.1 nG. A better understanding 
of the interplay between dynamo theory and the global conservation laws of magnetized plas- 
mas has been recently achieved [2] also because of the improved comprehension of the solar 
dynamo action. It is then plausible, within the dynamo hypothesis, that the protogalactic 
field could be even much smaller than the nG and still explain some crucial properties of our 
magnetized Universe. 

Thanks to magnetic flux (and magnetic helicity) conservation in a conductive plasma, 
a magnetic field of nG strength at the epoch of galaxy formation can be as large as mG 
(i.e. roughly 6 orders of magnitude larger) at the epoch of photon decoupling, i.e. for 
-2dcc = 1100. If large-scale magnetic fields have primordial origin, they were present prior to 
matter-radiation equality affecting, potentially, CMB anisotropics [1]. Through the years, 
various studies have been devoted to the effect of large-scale magnetic fields on the vector 
and tensor CMB anisotropics [3] (see also [4] and references therein for some recent review 
articles) . 

The implications of fully inhomogeneous magnetic fields on the scalar modes of the ge- 
ometry remain comparatively less explored. By fully inhomogeneous we mean stochastically 
distributed fields that do not break the spatial isotropy of the background [4]. One of the 
aims of the present paper is to partially bridge this gap and to open the way for further 
developments. In short the idea is the following. The simplest set of initial conditions for 
CMB anisotropics, implies, in a ACDM framework, that a nearly scale-invariant spectrum 
of adiabatic fluctuations is present after matter-radiation equality but before decoupling for 
typical wavelengths larger than the Hubble radius at the corresponding epoch [5]. It became 
relevant, through the years, to relax the assumption of exact adiabaticity and to scrutinize 
the implications of a more general mixture of adiabatic and non-adiabatic initial conditions 
(see [6] and references therein). In this paper it will be argued, along a similar perspective, 
that large-scale magnetic fields slightly modify the adiabatic paradigm so that their typical 
strengths may be constrained. To achieve such a goal, the first step is to solve the evolution 
equations of magnetized cosmological perturbations well before equality. The second step is 
to follow the solution through equality (and up to decoupling). On a more technical ground, 
the second step amounts to the calculation of the so-called transfer matrix whose specific 
form is one of the the subjects of the present analysis. 

Consider then the system of cosmological perturbations of a flat Friedmann-Robertson- 
Walker (FRW) Universe, characterized by a conformal time scale factor a{T), and consisting 
of a mixture of photons, baryons, CDM particles and massless neutrinos. In the conformally 
Newtonian gauge [6, 7] the scalar fluctuations of the metric tensor g^^, are parametrized in 
terms of the two longitudinal fluctuations i.e. 6goo = 2a^0(r, x) and 6gij = 2a^-?/'(r, x)(5jj. 
The Hamiltonian and momentum constraints, stemming from the (00) and (Oi) components 



of the perturbed Einstein equations are: 

V^^ - m{n(t> + ^') = A7iGa^[5p, + 5pb], 5pB(r, x) = f ^^,\ , (1) 

W\n<P + ^') = -AnGa^ip, + p,)e,, (2) 

where 7i = a' /a and the prime denotes a derivation with respect to the conformal time 
coordinate r. The total energy and pressure densities of the mixture, i.e. pt = I^aPa and 
Pt = SaPa, determine the evolution of the background geometry according to Friedmann 
equations: 

U^ = ^a^p^^ n^_n' = AnGa\p,+p,), p[ + 3n{pt + Pt) = . (3) 

In Eqs. (1) and (2) 6pt and 6^, denote, respectively, the total density fluctuation of the fluid 
mixture and the divergence of the total velocity fleld (i.e. 6t = divl) whose expressions, in 
terms of the four components of the plasma, i.e. z/, 7, c (CDM) and b (baryons), is 

^Pt = Yl '^Pa, Spt = Y^ 6p^, {pt + pt)Ot = ^(pa + Pa)^a- (4) 

The spatial components of the perturbed Einstein equations, imply, instead 

< + 7^(0' + 2^') + {2n' + re)(j) + \v\(t) - V^) = A'KGa'idp, + (5pb), (5) 

V^(0 - V^) = l2^Ga\p, + p.)VV, + {p^ + p^)VVb], dp^ = ^. (6) 

In Eq. (6) V^o-jy is the neutrino anisotropic stress, while V^ctb is the magnetic fleld anisotropic 
stress deflned as: 

where, VI-q{x) is the magnetic energy density referred to the photon energy density and it 
is constant to a very good approximation if magnetic flux is frozen into the plasma element 
[1, 3, 4], as assumed throughout the paper. The induced Ohmic current J is solenoidal (in 
the magnetohydrodynamical (MHD) description adopted here) and it is simply given by 
AtiJ = V X B. Moreover, in MHD, E = J/a ^ (V x B)/a (where a is the conductivity). 
Since, prior to decoupling, the Universe was a rather good conductor [4], the contribution 
of the electric energy density and of the Poynting vector appearing, in principle, in Eqs. 
(1) and (2) can be safely neglected for typical length scales much larger than the screening 
length of the plasma. It should be stressed that, in Eq. (6), on top of the magnetic piece, 
the only contribution to the anisotropic stress of the fluid mixture comes from massless 
neutrinos^ (that are collisionless for temperatures smaller than 1 MeV) and it is parametrized 

^If neutrinos would have a mass in the nicV range, they would be non-rclativistic today but they will 
still be counted as radiation prior to decoupling. 



by (Tj/. The evolution of 5px. can be determined from the covariant conservation of the (total) 
energy-momentum tensor: 

5p[ - 3^'(pt + Pt) + (Pt + Pt)^t + m(l + cl)5p, + 3Wp„ad = ^, (8) 

where c^ = p[/ p[ is the (total) sound speed and where the (total) pressure density fluctu- 
ation 5pt has been slpit into the adiabatic contribution (i.e. c^Spt) supplemented by the 
non-adiabatic pressure density fluctuation (i.e. (5pnad)- The electromagnetic contribution 
appearing in Eq. (8) contains an electric field and it is therefore suppressed. 

The evolution of the CDM component feels indirectly the presence of the magnetic field 
intensity through the Hamiltonian constraint (1) and the relevant equations are 

e'^ + n9c + V^(j) = 0, 5'^ = S^P' - Oc, 5^ = —. (9) 

Pc 

The neutrinos are coupled to the magnetic field through the Hamiltonian constraint (1) and 
through Eq. (6) (involving the neutrino anisotropic stress V^o";^): 

where, in full analogy with Eq. (9) the neutrino density contrast 5i, has been introduced. 
Photons and baryons are tightly coupled by Thompson scattering and form, effectively, a 
single fluid characterized by a velocity field ^^b = 0-y = 9\,. The relevant evolution equations 
are, in this case. 



0',. + TT^^^.^ + 77^^ + v^0 = 7 ,:;!::^::\ , (12) 



TiB^^ V% , 3 V-[JxB] 

;i + i?b) ^"^ 4(l + i?b) 4aV7(l + i?b 

where i?b is the baryon to photon ratio that depends on the redshift z. Deep in the radiation- 
dominated epoch i.e. for r <^ Tcq the solution for the magnetized adiabatic mode can be 
obtained, in Fourier space, by solving, simultaneously, Eqs. (9), (10) and (11)-(12). The 
compatibility of the magnetized adiabatic mode with Eqs. (l)-(2) and also with Eqs. (5)-(6) 
fixes the integration constants. Defining as k the Fourier (comoving) wave-number we shall 
be interested in wavelengths much larger than the Hubble radius, i.e. kr < 1. For r <^ Teq, 
the density contrasts for the magnetized adiabatic mode are, to lowest order in kr < 1, 

3 3 

5^ = 6,, = -2(j)i - R^^Ib, S^ = Sc = --(f)i- -R-y^B, (13) 

where the fractional contribution of photons to the radiation plasma, i.e. R^ has been 
introduced and it is related to R^, i.e. the fractional contribution of massless neutrinos, as 

R, = i-R., R^ = Y^^ ""=8^^(11) ^^■^^\-y)- ^^^^ 



From Eqs. (9), (10) and (12) the velocity fields of the various species are 









k'r^a^. (15) 



The quantities ipi{k) and (pi{k) appearing in Eqs. (13) and (15) denote the super-Hubble 
fluctuations that are initially present prior to equality. By solving in terms of the neutrino 
anisotropic stress a,, and by recalling Eq. (6) the relation between ipi and 0i can be obtained: 

^i = 0if 1 + lR^ + ^(4aB - R^Qb), ^u = -|^^B + ^(^i - 0i). (16) 

V 5 / 5 R^ bRu 

In the limit ctb ^ and VL-q -^ this solution reproduces the standard adiabatic mode 
in the longitudinal gauge (see third and fourth references in [7]). To follow the fate of the 
magnetized adiabatic mode through Teq it is practical to exploit the total density contrast 
on uniform curvature hypersurfaces (conventionally denoted by (^ [7, 8]) or the curvature 
perturbation on comoving orthogonal hypersurfaces (conventionally denoted by TZ [7, 8]) 
whose speciflc definitions, in terms of the variables of the longitudinal gauge, are 

( = -^-7^ , 11 = -i) — -—, C = TZ+ ■ -. 17 

p[ W - W l2'nGa\pt + Pt) 

The first and second relations in Eq. (17) are the definitions of C, and TZ in terms of the 
conformally Newtonian variables. The third relation in Eq. (17) can be obtained by substi- 
tuting the definitions of C, and TZ back into Eq. (1) and by recalling the background relations 
(3). In the limit when the relevant wavelengths are all larger than the Hubble radius at 
the corresponding epoch, i.e. kr < 1, the third relation in (17) implies that 7^ ~ (^. The 
evolution equation for ( can then be obtained by inserting the definition of ( into Eq. (8); 
the result is 

C = -^Sp^^, + -^ U - 1)SPB - |. (18) 

Pt + Pt Pt + Pt\ 3/ 3 

The non-adiabatic pressure density variation 6pna,d can be written as a sum of the relative 

entropy fiuctuations over the various components of the mixture 

5pnad = XTTT Yl Pi Pj(Csi " Csj)'^iJ' '^^ = ~3(Ci - Cj), C^i = 4' (19) 

O'T-Pt ij Pi 

where iSij are the relative fiuctuations in the entropy density that can be computed, from 
Eq. (19), in terms of the density contrasts of the individual fiuids, i.e. 

^ PUu + P'^Cy + P'cCc + P'hCh . . SpB ._. 

Pt 3(pt + Pt) 



In the case of adiabatic (magnetized) initial conditions it can be easily verified that Cc = Cb = 
Cu = C7 so that Spnad = 0. Deep in the radiation-dominated epoch, for r <^ Teq, c^ -^ 1/3 
and, from Eq. (18), C' = 0, so that 

C = Ci ^ 7^i, Ci = -^0i(i + ^R.) - ^(4aB - R.nB). (21) 

When the Universe becomes matter-dominated, after Teq, c^ ^ and the second term at 
the right hand side of Eq. (18) does contribute significantly at decoupling (recall that 
for /i^fimattcr = 0.134, Tdcc = 2.36 Teq). Consequently, from Eq. (18), recalling that c^ = 
4acq/[3(3a + 4acq)], we obtain 

Cf = G-T^f^^, fiBf = l^Bi. (22) 

4(3a + 4acq) 

The inclusion of one (or more) adiabatic modes changes the form of Eq. (18) and, con- 
sequently, the related solution (22). For instance, in the case of the CDM-radiation non- 
adiabatic mode the relevant terms arising in the sum (19) are Scy = So, = Si where Si is 
the (constant) fluctuation in the relative entropy density initially present (i.e. for r ^ Teq). 
If this is the case 6paad = c^PcSi and Eq. (18) can be easily solved. The transfer matrix for 
magnetized CMB anisotropies can then be written as 

Mss MsB I I ^5; I . (23) 





M 



BB 




In the case of a mixture of (magnetized) adiabatic and CDM-radiation modes, we flnd, for 

Mcc^l, M^s^~\, TWcB-^, Mss^l, A^5B^0, (24) 

and A^BB -^ 1- Equations (23) and (24) may be used, for instance, to obtain the magnetized 
curvature and entropy fluctuations at photon decoupling in terms of the same quantities 
evaluated for r ^ Teq. A full numerical analysis of the problem conflrms the analytical 
results summarized by Eqs. (23) and (24). The most general initial condition for CMB 
anisotropies will then be a combination of (correlated) fluctuations receiving contribution 
from 5pnad and from the fully inhomogeneous magnetic fleld. To illustrate this point, the 
form of the Sachs- Wolfe (SW) plateau in the sudden decoupling limit will now be discussed. 
To compute the SW contribution we need to solve the evolution equation of the monopole 
of the temperature fluctuations in the tight coupling limit, i.e. from Eqs. (11) and (12), 



^^ ' l + i?b ^ ■ 3 l + i?b "' ■ l + i?b" 3" " 3(l + i?b) 



<' + T^^r^< + -^-r^ = 4^" + 7^^^' - o^V - .,, , „, (^B - 4aB). (25) 



In the sudden decoupling approximation the visibihty function, i.e. /C(t) = K'{T)e~'^'^'^^ and 
the optical depth, i.e. e"''*-^-' are approximated, respectively, by 5(r — race) and by 6'(r — race)- 
The power spectra oi (, S and Qb are 



Vcik) 



Aci^ 



rir — l 



Vs{k)=As[^ 

Kr, 



Vn{k) = He)^l,[^^ 



2e 



(26) 



where Ar, As and Qbl are constants and 



He) 



\2e 



4(6-£)(27r 
e{3-2e)r^{e/2y 



n 



B L 



Pbl 

P7 ' 






jj=a\r)p,ir). (27) 



To deduce Eqs. (26) and (27) the magnetic field has been regularized, according to a common 
practice [3, 4], over a typical comoving scale L = 27i/kL with a Gaussian window function 
and it has been assumed that the magnetic field intensity is stochastically distributed as 

{B,{k,T)B^{p,T)) = '^pi{k)PB{k,T)6^'\k + p), (28) 



k^ 



where 



Phk) 



S' 



ki^ 



PB{k,r) = A^ 



(29) 



As a consequence of Eq. (28) the magnetic field does not break the spatial isotropy of the 
background geometry. The quantity kp appearing in Eqs. (26) and (29) is conventional 
pivot scale that is 0.05Mpc(see [6] for a discussion of other possible choices). Equations 
(26) and (27) hold for < e < 1. In this limit the magnetic energy spectrum is nearly 
scale- invariant. This means that the effect of the magnetic and thermal diffusivity scales 
(related, respectively, to the finite value of the conductivity and of the thermal diffusivity 
coefficient) do not affect the spectrum [4]. In the opposite limit, i.e. e ^ 1 the value of the 
mode-coupling integral appearing in the two-point function of the magnetic energy density 
(and of the magnetic anisotropic stress) is dominated by ultra-violet effects related to the 
mentioned diffusivity scales [4]. Using then Eqs. (26) and (27) the Ce can be computed for 
the region of the SW plateau (i.e. for multipoles i < 30): 



G 



25 ^ 



9 
100 



— — -R fiBL^2(e, 



- —^jAi^AsZi{nrs,i) cos7rs 



4 3 / — — 

+—As2:i{ns,i) - —^JA^RyfiBL23{nr,e,i)cos%r 

6 / — — 
+— y As R-^/^BL23{ns,e,i) cos -fbs 



(30) 



where the functions Zi, Z2 and Z^, 



Zi{n,- 



4 U. 



„-i r(3 

2"- 



n)r{i + 



n—1 



r2(2- f ir 



_|_ o n 

'2 2 



(31) 



2 VA;^/ r2(|-£)r(£ + 2-£) 



^3(n,M) ^ ^2-2-^v^rferrfer 'r''Tr^^^'\ ,(33) 



4 VA:^y VA:py r2fl-|-Mr(! + 



£ 

4 2 4/"V4''' 2 

are defined in terms of tlie magnetic tilt e and of a generic spectral index n which may 
correspond, depending on the specific contribution, either to n^ (adiabatic spectral index), 
or to ns(non-adiabatic spectral index) or even to rirs = {ur + ns)/2 (spectral index of the 
cross-correlation). In Eq. (30) 7rs, 7fer and 7^6 are the correlation angles. In the absence 
of magnetic and non-adiabatic contributions and for Eqs. (30) and Eq. (31) imply that for 
n-r = 1 (Harrison-Zeldovich spectrum) i{i + l)Q/27r = A(;/25 and WMAP data [5] would 
imply that A,^ = 2.65 x 10~^. Consider then the physical situation where on top of the 
adiabatic mode there is a magnetic contribution. If there is no correlation between the 
magnetized contribution and the adiabatic contribution, i.e. •jbr = 7r/2, the SW plateau 
will be enhanced in comparison with the case when magnetic fields are absent. The same 
situation arises when the two components are anti-correlated (i.e. cos7br < 0). However, 
if the fluctuations are positively correlated (i.e. cos '-fbr > 0) the cross-correlation adds 
negatively to the sum of the two autocorrelations of ( and Q^ so that the total result may 
be an overall reduction of the power with respect to the case 75^ = vr/2. In Eq. (31), (32) 
and (33) ko = Tq^ where tq is the present observation time. Taking as typical parameters 
h^n^ = 0.111, h'^Qk = 0.023 and h'^Q, = 4.15 x lO'^ (with h = 0.73), the amplitude of the 
magnetic field intensity can be constrained by requiring that the adiabatic mode dominates. 
In the case when the magnetic and adiabatic contribution are totally anticorrelated (i.e. 
cos '-fhr = ~1)) which is the most unfavourable case the bound on the protogalactic field 
read, at the present epoch and over a typical comoving scale L = 271 /k^ with ki = Mpc~^ 

Bl < 2.5 X 10"^ G, rir ~ 0.951, e ~ 0.9. (34) 

As indicated, the bound (34) assumes a nearly scale- invariant (but slightly red [5]) adiabatic 
mode and the maximally allowed magnetic spectral tilt. A further reduction of e leads to a 
slight relaxation of the bound; for instance for e ^ 0.4, B^ < 6.3 x 10~^ G. 

In CMB physics is common practice to perform model-independent analysis on the pa- 
rameter space of the allowed initial conditions by including, for instance, correlated (or an- 
ticorrelated) non-adiabatic modes in the game (see, for instance, [6] and references therein). 
Up to now the effects related to fully inhomogeneous magnetic fields have been discussed 
within a different standard whose limitations were the impossibility of defining accurately 
both initial conditions and normalization of the magnetized (scalar) CMB anisotropics. The 
results reported here allow to overcome these difficulties and lead naturally to a strategy of 
parameter extraction where large-scale magnetic fields are treated consistently as a further 
degree of freedom in the space of the initial conditions. According to this perspective, it 

7 



will be important to pursue the analysis of small-scale effects by semi-analytical methods to 
corroborate and interpret more numerical studies related to parameter extraction. The first 
step in this direction would be to go to higher orders in the tight coupling expansion and 
generalize the standard treatment to the case when the pre-decoupling plasma is effectively 
magnetized. Along these directions work is in progress. 
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